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Timing of inorganic phosphate release modulates 
the catalytic activity of ATP-driven rotary motor 
protein 

Rikiya Watanabe 1,2 & Hiroyuki Noji 1 



F r ATPase is a rotary motor protein driven by ATP hydrolysis. The rotary motion of F r ATPase 
is tightly coupled to catalysis, in which the catalytic sites strictly obey the reaction sequences 
at the resolution of elementary reaction steps. This fine coordination of the reaction scheme 
is thought to be important to achieve extremely high chemomechanical coupling efficiency 
and reversibility, which is the prominent feature of F r ATPase among molecular motor 
proteins. In this study, we intentionally change the reaction scheme by using single-molecule 
manipulation, and we examine the resulting effect on the rotary motion of F r ATPase. When 
the sequence of the products released, that is, ADP and inorganic phosphate, is switched, we 
find that F-| frequently stops rotating for a long time, which corresponds to inactivation of 
catalysis. This inactive state presents MgADP inhibition, and thus, we find that an improper 
reaction sequence of F r ATPase catalysis induces MgADP inhibition. 
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F!-ATPase (a 3 p 3 y5s), a catalytic sub-complex of FoFi-ATP 
synthase, is a rotary motor protein fuelled by ATP 
hydrolysis 1-3 . Three subunits (oc 3 p 3 y) function as the 
minimum component of the rotating system, where the a 3 |3 3 
unit forms a cylindrical stator and the y rotor subunit penetrates 
the centre of the cylinder 4 . The catalytic sites for ATP hydrolysis 
are located on each interface between the a and P subunits, 
mainly on the P subunits, and the three P subunits are always in 
different catalytic states 4 . The interconversion of the catalytic 
states of P subunits, which induces their conformational 
transition, drives the rotation of the y subunit 5-7 . 

The rotary motion of F : can be directly visualized by optical 
microscopy 8-10 . The unitary step size of the rotation of F : is 120°, 
and each step is coupled to a single turnover of ATP hydrolysis 11 . 
The 120° step is further divided into 80° and 40° sub-steps 12 ' 13 . 
The 80° sub-step is triggered by ATP binding and ADP release, 
each of which occurs on different P subunits 14 . The 40° sub-step 
is triggered by ATP hydrolysis and release of inorganic phosphate 
(Pi), which also occurs on different P subunits 14 ' 5 . The catalytic 
state of the P subunits is tightly coupled to the rotation of the 
y subunit; thus, F : can synthesize ATP from ADP and P A when 
the y subunit is forcibly rotated in the reverse direction 16 ' 17 . 
Mechanically induced ATP synthesis is the physiological role of 
F : of FoFi-ATP synthase, in which F 0 , a proton-driven motor 
protein, compels F : to rotate in the reverse direction for ATP 
synthesis. This reversibility of chemomechanical coupling is the 
prominent feature of F : among other motor proteins. 

To achieve reversibility of chemomechanical coupling, fine 
coordination of the reaction sequence, which properly induces the 
conformational transition of P subunits, is required. To under- 
stand the fine coordination mechanism, extensive studies have 
been dedicated to unveiling the basic reaction scheme for the 
rotation and catalysis of F : , focusing on the catalytic turnover of 
individual P subunits (Fig. la). In this scheme, an individual P 
subunit binds to ATP when the y subunit is oriented at a specific 
angle, and the binding angles for the individual P subunits differ 
by ± 120°. Each P subunit hydrolyses the bound ATP into ADP 
and Pi after y rotates another 200° from the ATP-binding angle 18 . 
In addition, the produced ADP and Pi are released from 
the catalytic site after additional 40° and 120° rotations, 
respectively 14 ' 15 ' 19 ' 20 . When the y subunit returns to the 
original angular position, the P subunit initiates the next round 
of catalysis by binding to a new ATP. Thus, the reaction 
sequences among the three P subunits are well coordinated at the 
resolution of elementary reaction steps, which contributes to the 
reversibility of the chemomechanical coupling on F^ 

Then, the question arises as to what happens when the reaction 
sequence of F : is changed. In the case of the hydrolysis step, 
although the timing of hydrolysis shifts by 80° relative to the 
authentic reaction angle, F x drives the rotation without interrup- 
tion, demonstrating the robustness of the chemomechanical 
coupling of F : (ref. 21). In contrast, in the case of other 
elementary reaction steps, such as Pi release, the effect of timing 
on rotary motion has not been examined so far, although we 
established a method to change the timing of Pi release by using 
single-molecule manipulation 15 . In particular, the free energy 
change resulting from Pi release was relatively high among 
reaction steps; that is, the rotary torque was mainly generated by 
Pi release, and therefore, the effect of the timing of Pi release on 
rotary motion is expected to be larger than that of other reaction 
steps 4,1 5 ' 22 . In addition, the timing of Pi release in F : is different 
from that in other motor proteins fuelled by ATP hydrolysis, such 
as kinesin and myosin, in which P t release occurs before ADP 
release 23 ' 24 , and this difference is thought to be important for the 
reversibility of the chemomechanical coupling mechanism of 
Fi (ref. 15). 
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Figure 1 | Chemomechanical coupling scheme of FrATPase. (a) The 

circles and red arrows represent the catalytic state of the (3 subunits and the 
angular positions of the y subunit. Each p subunit completes one turnover 
of ATP hydrolysis in a turn of the y subunit, where the three P subunits vary 
in their catalytic phase by 120°. Regarding the catalytic state of the top p 
subunit (cyan), ATP binding, hydrolysis, ADP release and inorganic 
phosphate (Pi) release occur at 0°, 200°, 240° and 320°, respectively, 
(b) Time course of the rotation of the hybrid F r ATPase, a 3 p(E190D) p 2 y, 
in the presence of 1 mM ATP. Red, orange, and green represent the pause 
at 200°, 320°, and 360° for p(E190D), respectively. The inset shows a 
magnification of the time course. 

In this study, to further investigate the role of the finely 
coordinated reaction sequence of F^ we evaluate the effect of 
changing the timing of Pi release on the chemomechanical 
coupling mechanism. When is released before ADP, we find 
that F : frequently stops rotating for a long time, which 
corresponds to inactivation of catalysis. This inactive state 
presents MgADP inhibition, and thus, we find that an improper 
timing of Pi release induces MgADP inhibition. 

Results 

Rotation of hybrid ¥ v We observed the rotation of hybrid F : 
carrying one |3(E190D), that is, a 3 p 2 |3(E190D)y, in the presence 
of 1 mM ATP. Glul90 of the |3 subunit, the so-called 'general 
base', is known to be one of the most important residues for 
hydrolysis 25-27 , and its substitution into aspartic acid causes 
distinctively slow hydrolysis of ATP 28 and strong temperature 
sensitivity . Therefore, the kinetic steps of the mutated catalytic 
site are easily distinguished from those of the other two sites. In 
the presence of 1 mM ATP, the hybrid F : showed three distinctive 
dwells caused by the incorporated |3(E190D); the dwell for the 
temperature- sensitive reaction at 0° (t = 32 ms) 15 ' 29 , which is the 
same as the dwell at the ATP-binding angle, ATP hydrolysis 
at 200° (t = 318 ms), and P t release at 320° (i = 9ms) (cyan in 
Fig. la,b). The dwell time for P^ release at 320° was prolonged 
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Figure 2 | Single-molecule manipulation for Pj release, (a) Schematic representation of manipulation procedures. When the hybrid F r ATPase paused at 
200°; that is, the angle for ATP hydrolysis by p(E190D), we switched on the magnetic tweezers and arrested Ft at ±10° from 200°. After release from 
arrest, F r ATPase roughly showed two behaviours: rotating forward to next reaction angle CON'), or staying at the original pausing angle ('OFF'). In 
addition, the behaviour of 'ON' was classified into three types: (I) rotating forward to 320° and showing the distinctive pause for P, release (time constant 
(t) is ~9ms); (II) rotating forward to 360° (skip the pause at 320°); or (III) rotating forward to 320° and showing a long pause (t= ~39s). 
(b) Hydrolysis-waiting time at 200° before F r ATPase is arrested with the magnetic tweezers. The average waiting time was 370 ± 37 ms for the 
behaviours (I) and (II) (light grey), and 275 ± 56 ms for the behaviour (III) (grey). Error bars indicate s.d. (c) Time course of P Q n- Red and blue points 
represent P Q n in the absence or presence of 100 mM Pj, respectively. The number of arrest trials and molecules used for each data point were 23-108 
and 8-19. The time courses were fitted with the reaction scheme; F-|-ATP<-> -F-| (ADP + P i )->F 1 -ADP + P, (solid lines), (d) Example of the time course 
of the arrest experiments (blue periods) corresponding to the behaviour of (I), (II) and (III), respectively, (e) Magnification of d. Orange and green 
represent the pauses at 320° and 360°. 



because of the viscous drag of the rotary probe on y as previously 
reported 30 ' 31 . Hereafter, we refer to the angle for ATP hydrolysis 
by P(E190D) as the catalytic angle. 



Single-molecule manipulation to induce Pi release. During free 
rotation, one |3 hydrolyses ATP and another |3 releases the pro- 
duced Pi at the catalytic angle, resulting in a rotation of the y 
subunit in an anticlockwise direction (Fig. la). In contrast to the 
free rotation, we recently succeeded in inducing the release of 
another produced Pi from the same catalytic site used for ATP 
hydrolysis by arresting the y subunit for a long time ( > 10 s) 15 ' 22 . 



In this study, we used the same method to induce the release of 
produced Pi immediately after ATP hydrolysis. For manipulation 
of y rotation, a magnetic bead was attached to the y subunit of F l5 
and the oc 3 p 3 ring was immobilized on the glass surface. When the 
hybrid F : showed a pause for ATP hydrolysis by |3(E190D) at 
200°, we turned on the magnetic tweezers within ~ 450 ms 
(Fig. 2a,b) and arrested F : at ± 10° from 200°. After a certain 
time period had elapsed, we turned off the magnetic tweezers and 
released F : from arrest. Then, F : roughly showed one of two 
behaviours without exception, as previously reported: rotating 
forward or restarting the hydrolysis-waiting pause (Fig. 2a). 
When F : shows the former behaviour, |3(E190D) has already 
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Figure 3 | Analysis of the effect of Pj release at 200°. (a,b) Histograms of the dwell time at 320° after the release from 200° arrest for 10 s, 30 s, 

and 60s, or during free rotation. The bin sizes of the histogram are 15s for a and 4ms for b. Blue bars in a represent the number of pauses longer 
than 100 ms (behaviour (III)). The data after first bin in b were well-fitted by the single exponential curve where i = 9ms (P, release dwell at 320°). 
The deviations of the first bins from the single exponential curves are coloured green (behaviour (II)). (c) Probability of Pj release at 200°. The red solid 
line presents the simulation curve for P, release calculated from the analysis of P Q n based on the kinetic scheme as follows: F r ATP<-> F-i • (ADP + Pj)-> 
F r ADP + Pj. Green, blue and orange points represent the probability of occurrence of II (P M ), III (P m ) and their sum (Pn + Pm), respectively. P M and P m 
were determined from the analysis of b and a, respectively. 
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Figure 4 | Analysis of the long pause at 320°. (a) Time course of P m in the absence (red) or presence of 100 mM Pj (blue). The red solid line is the 
fitting curve based on the kinetic scheme [1]. The blue solid line is the simulation curve based on the kinetic parameters determined from the fitting of 
Pm in the absence of Pj. (b) Relationship between P m for 60s arrest and Pj concentration. The simulation curve in Fig. 4b was redrawn against the 
P, concentration, (c) Histogram of durations of the long pauses at 320° from the experiments of 10-60 s arrests. Solid line represents the fitting with 
single exponential curve where the time constant (Tj_> a ) is 39 s. Bin size is 20 s. 



hydrolysed ATP and exerted a torque on the magnetic beads, 
whereas when P(E190D) shows the latter behaviour, it has not 
hydrolysed ATP because F : cannot generate torque unless 
P(E190D) catalyses hydrolysis. These behaviours are hereafter 
referred to as 'on' and 'off, respectively. Then, we measured the 
probability of the post-hydrolysis state as the probability of the 
'on' case against the total number of trials, Pqn- In the absence of 
Pi in solution, Pqn increased toward 100% as the arrest time 
increased; in contrast, in the presence of 100 mM Pj, Pqn did not 
change depending on the arrest time, which suggests that Pqn 
depends on the time course of Pj release from the catalytic site 
(Fig. 2c). 

We extensively analysed the 'on' case and found that it can be 
classified into three behaviours; (I) rotating forward to 320°, 
(II) rotating forward to 360°; that is, skipping the pause at 320°, 
and (III) rotating forward to 320° and spontaneously stopping the 
rotation for a long time (Fig. 2d,e). When |3(E190D) shows (I), 

4 



it has hydrolysed ATP during its arrest but has not released 
the produced Pi because the rate-limiting step at 320° is not 
hydrolysis by (3(WT) but rather P 4 release by P(E190D) (Fig. la). 
F : cannot generate torque at 320° unless P(E190D) releases Pi. 
When P(E190D) shows (II), it has hydrolysed ATP and 
simultaneously released the produced Pi, and therefore, F : has 
already exerted torque for the step from 320° to 360° driven by P t 
release. Regarding the behaviour of (III), we have not identified 
the chemical state of P(E190D) thus far; however, the pausing 
angle is the same as that of the MgADP- inhibited state . 
Hereafter, we characterize (III). First, we analysed the dwell time 
at 320° and calculated the probability of occurrence for each 
behaviour. As shown by the blue bars in Fig. 3a, distinctively long 
dwell times (longer than 100 ms), which correspond to the 
behaviour of (III), can be recognized. Then, we calculated the 
probability of occurrence of (III), that is, P m , as the number of 
events corresponding to the blue bars divided by total number of 
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Figure 5 | Mechanism of MgADP-inhibited state. Black solid line: the 
normal reaction pathway of active F r ATPase. FrATPase releases P, at 320°. 
Grey lines: alternative reaction pathways. F r ATPase releases Pi before ADP 
at 200°. F r ATPase lapses into MgADP inhibition at 320° with some 
probability (dashed line). 

events in the histograms, and found that P m increased as the 
arrest time increased. As shown in Fig. 3b, we changed the time 
scale of the histogram to shorter bins (4 ms), and found that the 
histogram deviated from a single exponential decay model as the 
arrest time increased. For example, at 60 s of arrest, the first bin 
(<4ms) contained many more events in which |3(E190D) 
skipped the pause at 320° than expected from a single 
exponential decay that fit the data well after the first bin with a 
time constant of 9 ms (black line in Fig. 3b), which was consistent 
with that of Pi release in free rotation (rightmost panel in Fig. 3b). 
Therefore, the first bin deviated from single exponential decay 
because of the behaviour of (II); accordingly, the first bin 
potentially includes behaviours of both (I) and (II). To calculate 
the probability of occurrence of (II), that is, P u , we fit the data 
points of the histograms (except for the first bin) to a single 
exponential curve with 9 ms as the time constant of Pi release. 
P n was determined as the number of events that deviated from 
the single exponential curve in the first bin (green bars in Fig. 3b) 
divided by total number of events of the histogram in Fig. 3 a. In 
addition, we also calculated the probability of occurrence of (I), 
that is, Pi, as the number of remaining events (light yellow bars in 
Fig. 3b) divided by total number of events. 

Kinetic analysis of Pi release and the inactive state. From the 
analysis of Pqn based on the irreversible reaction scheme (Fig. 2c) 
F : • ATP Fx • (ADP + Pi) -> Fi • ADP + Pi, we can calculate the 
probability of Pi release during arrest as P^ff, as shown by the 
red line in Fig. 3c. For example, P^ff is 37% and 60% for 30-s and 



60 -s arrest, respectively. Here, we compared P n and P^ff, an d 
found that P n was much smaller than Poff (Fig- 3c, green). 
As mentioned above, we could not identify the chemical state of 
(III). Thus, we suspected that (III) was caused by Pi release 
immediately after ATP hydrolysis and summed P n and P m . The 
sum of Pu and P m was almost the same as P^ff (Fig. 3c, orange), 
suggesting that (III) was caused by Pi release immediately after 
ATP hydrolysis. To confirm the relationship between (III) and 
Pi release, we conducted the arrest measurement in the presence 
of high [Pi] (~100mM), where Pi did not affect the kinetics of 
ATP binding and produced ADP release 14,33 . When we added 
more than lOmM Pi in solution, P m was completely suppressed 
toward 0% (Fig. 4a,b), suggesting that (III) was induced after the 
Pi release. Then, the kinetic scheme of (III) is given below 
(equation (1)), where, k^yd* ^I P > ^off> ^oA> an d are the rate 
constants of ATP hydrolysis, synthesis, Pi release, binding and 
inactivation of P(E190D), respectively. 

k ATP k Pi 

Fi • ATP ^± Fi • (ADP + Pi) Fi • ADP + P t 

Kyn *S 

^F!-ADP* + Pi (1) 

As previously reported, /chyd P > ^£ P > ^off> an d fcon were given as 
2.3 s _1 , 0.8 s _1 , 0.021 s _1 and 2 x 10 3 M~ 1 s~ \ respectively 15 . 
Then, from the fitting of P m in Fig. 4a based on the kinetic 
scheme as shown above, the only fitting parameter; fc a _>i, was 
determined as 0.019 s _1 . To verify this kinetic scheme, we 
fixed the arrest time as 60s, and we re-plotted P m against P t 
concentration. Based on this kinetic scheme with the determined 
/c a ^i, we mathematically reproduced the experimental result of 
P m as shown in Fig. 4b. From these results, we confirmed that 
this inactive state was induced after Pi release at 200° with a rate 
constant of 0.019 s _1 . 

We analysed the dwell time for activation from this inactive 
state (Fig. 4c). As shown in Fig. 4c, the histogram of the dwell 
time fit well to a single exponential curve, and we determined the 
time constant of activation, Ti^ a , as 39 s, which was almost the 
same as that for activation from an MgADP-inhibited state 32 . In 
addition, the aforementioned phenomena in which the addition 
of a large amount of P t prohibited F : from lapsing into the 
inactive state (Fig. 4a) is a well-known feature of the MgADP- 
inhibited state 34 , suggesting that the inactive state revealed in this 
study corresponds to MgADP inhibition. 

Discussion 

MgADP inhibition, a common feature of FrATPase and ATP 
synthases from various sources 32 ' 34-37 , is the catalytically inactive 
state. Although many biochemical studies have addressed 
MgADP inhibition, the mechanism of MgADP inhibition has 
remained elusive in the context of catalysis -rotation scheme. 
Here, we propose a simple model for how F x lapses into the 
MgADP-inhibited state during catalysis (Fig. 5). We assume that 
the fundamental principle of MgADP inhibition is the loss of 
the driving energy of rotation derived from Pi release. Our 
experimental data indicate that F : can continue rotating even if it 
releases Pi before ADP at 200°, suggesting that F : principally 
stores the driving energy derived from Pi release in conformation 
for 320°-360° rotation (grey solid line in Fig. 5). However, with a 
time constant of 0.019 s ~ l , F : dissipates the stored energy during 
rotation before reaching 320° and lapses into MgADP inhibition 
at 320° (grey dash line in Fig. 5). The driving force of 320-360° 
rotation after inhibition is a thermal agitation, as our previous 
work suggests that thermally agitated rotation of y over 30° 
induced activation from the MgADP-inhibited state 37 . Thus, the 
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MgADP- inhibited state is so stable that it takes a long time to 
resume rotation, Ti^ a = 39s (Fig. 4c) 32 . 

The crystal structures co- crystallized with azide, including the 
original structure 4 ' 38 , revealed that P at 320° (|3 E ) does not bind 
nucleotides or phosphate. Considering that azide is well known 
to stabilize the MgADP -inhibited state of these structures 
demonstrate the structural features of the MgADP-inhibited state. 
Thus, it is highly probable that P in the 320° state of MgADP- 
inhibited F : has no bound nucleotide or phosphate, which is 
consistent with our model (Fig. 5). 

It should be mentioned that the Walker group has challenged 
our reaction scheme of catalysis and rotation for active F ls on 
which the present work is based, that is, where Pi is released from P 
in 320° state. Based on the recent crystal structure of MP 1 in which 
P in the 320° state is bound to Mg-free ADP 39 , they proposed that 
Mg and P A are released before ADP. Although the reason for this 
apparent discrepancy is not clear, it should be noted that most of 
the crystal structures solved so far had no bound nucleotide on P in 
the 320° state; the only exceptions are the above-mentioned 
structure 39 and a structure in a different study 40 . In addition, there 
are several experimental differences between the present study and 
the crystal structure study such as the species from which F : is 
derived and the buffer contents. A precise correlation between the 
crystal structure and the reaction model based on single-molecule 
experiments remains to be established. 

In vivo, ¥i forms FoFi-ATP synthase and synthesizes ATP by 
coupling with a rotary motion driven by F 0 . Previous studies 
revealed that MgADP inhibition suppresses ATP hydrolysis but not 
synthesis 41 . Recently, the rotary motion of the ¥ 0 ¥i complex coupled 
to the synthesis of ATP was visualized at the single-molecule level 42- 
44 For further understanding of the effect of MgADP inhibition on 
ATP synthesis, verification by single-molecule measurement of the 
¥ 0 F 1 complex is desirable in the future. 

Methods 

Rotation assay. To visualize the rotation of hybrid F l5 a 3 (3(E190D) (3 2 y, the stator 
region (oc 3 (3(E190D)(3 2 ) was fixed on the glass surface, and a magnetic bead 
(4> = 0.3 urn; Seradyn, USA) was attached to the rotor part (y) as a probe for 
rotation and manipulation. The experimental procedure was as follows. The flow 
chamber was constructed from an uncoated top coverslip and a bottom coverslip 
whose surface was modified with Ni 2 + NTA. F x solution, which was diluted in 
buffer A (50 mM MOPS-KOH, 50 mM KC1 and 5mM MgCl 2 , pH 7.0) to a final 
concentration of 200 pM, was infused into the flow chamber. After 5 min, unbound 
F x molecules were washed out with buffer A containing 10% bovine serum 
albumin, and then streptavidin- coated magnetic beads in buffer A were infused. 
After 10 min, unbound beads were washed out with buffer A. Finally, buffer A was 
infused with the indicated amount of ATP or Pi. The rotating beads were observed 
under a phase-contrast microscope (IX-70 or IX-71; Olympus, Japan) with a x 100 
objective lens. The rotation assay was performed at 23 °C. 

Manipulation with magnetic tweezers. The stage of the microscope was 
equipped with magnetic tweezers that were controlled using custom-made software 
(Celery, Library, Japan). The rotary motion of the bead was recorded at 30 and 
3,000 frames per second, simultaneously (FC300M, Takex, Japan; FASTCAM 
1024PCI-SE, Photron, Japan). Images were stored in the hard disk drive of a 
computer as an AVI file and analysed using custom-made software. From the 
recorded images at 3,000 frames per second, we judged whether the operations 
precisely arrested the y subunit at the targeted angle or not. 
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